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The forward electroproduction of two light vector mesons is the first 
example of a collision process between strongly interacting colorless parti- 
cles for which the amplitude can be written completely within perturbative 
QCD in the Regge limit with next-to-leading accuracy. This amplitude can 
be written as a convolution of two impact factors for the virtual photon 
to light vector meson transition with the BFKL Green's function. In this 
lecture we first describe how the relevant impact factor is calculated, then 
we perform the convolution with the BFKL Green's function and illustrate 
the numerical procedure to obtain a well-behaved amplitude. 

PACS numbers: 12.38.Bx, 13.60.Le, 11.55.Jy 



1. Introduction 

In the BFKL approach [I], both in the leading logarithmic approxima- 
tion (LLA), which means resummation of leading energy logarithms, all 
terms (a s ln(s)) n , and in the next-to-leading approximation (NLA), which 
means resummation of all terms a s (a s ln(s)) n , the (imaginary part of the) 
amplitude for a large-s hard collision process can be written as the convo- 
lution of the Green's function of two interacting Reggeized gluons with the 
impact factors of the colliding particles (see, for example, Fig. [1]). 

* Presented by A. Papa at the "School on QCD, Low- a; Physics, Saturation and Diffrac- 
tion", Copanello (Calabria, Italy), July 1 - 14, 2007. 
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Fig. 1. Schematic representation of the amplitude for the 7*(p)7*(p') 
V(pi) V(p2) forward scattering. 



The Green's function is determined through the BFKL equation. The 
NLA singlet kernel of the BFKL equation has been achieved in the forward 
case [2] , after the long program of calculation of the NLA corrections [3] (for 
a review, see Ref. [1]). For the non- forward case the ingredients to the NLA 
BFKL kernel are known since a few years for the color octet representation 
in the t-channel [5]. This color representation is very important for the check 
of consistency of the s-channel unitarity with the gluon Reggeization, i.e. 
for the "bootstrap" [HE]- Recently it was completed also the calculation of 
the non-forward NLA BFKL kernel in the singlet color representation, i.e. 
in the Pomeron channel, relevant for physical applications [8]. 

On the other side, NLA impact factors have been calculated for colliding 
partons [91 [10] and for forward jet production [llj. Among the impact 
factors for transitions between colorless objects, the most important one 
from the phenomenological point of view is certainly the impact factor for 
the virtual photon to virtual photon transition, i.e. the 7* — > 7* impact 
factor, since it would open the way to predictions of the 7*7* total cross 
section. Its calculation is rather complicated and it was completed after 
year-long efforts [T2l [T3l [Ti] . 

A considerable simplification can be gained if one considers instead the 
impact factor for the transition from a virtual photon 7* to a light neutral 
vector meson V = p ,uj,cj). In this case, indeed, a close analytical expres- 
sion can be achieved in the NLA, up to contributions suppressed as inverse 
powers of the photon virtuality [15]. In particular, it turns out that (a) the 
dominant helicity amplitude is that for the transition from longitudinally 
polarized virtual photon to longitudinally polarized vector meson; (b) the 
impact factor, both in the LLA and in the NLA, factorizes into the convo- 
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lution of a hard scattering amplitude, calculable in perturbative QCD, and 
a meson twist- 2 distribution amplitude |15j . 

The knowledge of the 7* —* V impact factor allows for the first time 
to determine completely within perturbative QCD and with NLA accuracy 
the amplitude of a physical process, the 7*7* — ► VV reaction [I6| IT7]. This 
possibility is interesting first of all for theoretical reasons, since it can shed 
light on the role and the optimal choice of the energy scales entering the 
BFKL approach. Moreover, it can be used as a test-ground for comparisons 
with approaches different from BFKL, such as DGLAP, and with possible 
next-to-leading order extensions of phenomenological models, such as color 
dipole and /ct-factorization. But it could be interesting also for the possible 
applications to the phenomenology. Indeed, the calculation of the 7* — > V 
impact factor is the first step towards the application of BFKL approach 
to the description of processes such as the vector meson electroproduction 
7*p — * Vp, being carried out at the HERA collider, and the production of 
two mesons in the photon collision, 7*7* — > VV or 7*7 — > VJ/^f, which 
can be studied at high-energy e + e~ and colliders. 

In this paper we concentrate on the NLA forward amplitude for the 
7*7* — > VV reaction (Section 2). Such a process has been studied in 
Ref. |18] in the Born (2-gluon exchange) limit for arbitrary transverse mo- 
mentum and, for the forward case only, in Ref. [19] with LLA plus an 
estimate of NLA effects!]] 

First of all, we show how the available results for the 7* — > V impact 
factor (Section 3) and the BFKL Green's function can be put together to 
build up the NLA amplitude of the 7*7* — * VV process in the MS scheme 
(Section 4). Then we restrict ourselves to the particular case of collision 
of virtual photons with equal virtualities and present some numerical es- 
timates of our result, aimed at showing the extent of the contributions to 
the NLA amplitude from the impact factor and from the NLA kernel and 
the dependence on the energy scale introduced in the BFKL approach and 
on the renormalization scale which appears in the MS scheme. We show 
that, despite being the NLA corrections large and of opposite sign with re- 
spect to the leading order, it is possible to achieve a well-behaved form of 
the amplitude, by a suitable choice of the energy and renormalization scale 
parameters (Section 5). 

Then, we compare different procedures to optimize the perturbative re- 
sult and different representations of the amplitude, in order to have an 
estimate of the systematic effects which underlie our determination (Sec- 
tion 6). Finally, we calculate the differential cross section at the minimum 
squared momentum transfer and compare it with the approach of Ref. |19j 



1 The QCD factorization properties of this amplitude have been studied In Ref. |20| . 
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(Section 7). 

The use in our approach of the BFKL kernel improved by the inclusion 
of subleading terms generated by renormalization group analysis, which 
has been suggested to cure the instabilities in the behavior of the BFKL 
Green's function in the next-to-leading approximation [21], has been studied 
in Ref. [22] and is presented in Ref. [23]. The use of such an improvement 
has allowed for the numerical determination of the amplitude also in the 
case of colliding photons with strongly ordered virtuality. 



2. The amplitude for the electroproduction of two light vector 
mesons: kinematics and BFKL structure 

We consider the production of two light vector mesons (V = p°, u, (j)) in 
the collision of two virtual photons, 

l*{p)l*{p')^V{ Pl )V{p 2 ) . (2.1) 

Here, p\ and p 2 are taken as Sudakov vectors satisfying p\ = p\ = and 
2(pip2) = s; the virtual photon momenta are instead 

p = api p 2 , p = ap 2 — pi , (2.2) 

as as 

so that the photon virtualities turn to be p 2 = —Q 2 and (p') 2 = —Q 2 . We 
consider the kinematics when 

s » Q\ a » A 2 QCD , (2.3) 

and 

a= l + Ql + 0(s- 2 ), a' = l + 9l + (s- 2 ). (2.4) 
s s 

In this case vector mesons are produced by longitudinally polarized pho- 
tons in the longitudinally polarized state [15j . Other helicity amplitudes 
are power suppressed, with a suppression factor ~ mv/Qi,2- We will dis- 
cuss here the amplitude of the forward scattering, i.e. when the transverse 
momenta of produced V mesons are zero or when the variable t = (p\ —p) 2 
takes its maximal value to = — QiQi/ s + C( s_2 )- 

The forward amplitude in the BFKL approach may be presented as 
follows 

lm s (A) = / -^(gi,^) / $ 2 (-g 2)S0 ) (2.5) 

{Zir) J q 1 J q 2 



5+ioo 

X 
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This representation for the amplitude is valid with NLA accuracy. Here 
$i(q\, so) and ^ 2 (— Q2, sq) are the impact factors describing the transitions 
7*(p) — * V(pi) an d J*(p') — ► V(p2), respectively. The Green's function in 
(|2.5p obeys the BFKL equation 

$ 2 (qi ~ 02) = uG w {quq 2 ) - J d 2 qK(q 1 , q) G u (q, q 2 ) , (2.6) 

where K(q\,q2) is the BFKL kernel. The scale so is artificial. It is intro- 
duced in the BFKL approach at the time to perform the Mellin transform 
from the s-space to the complex angular momentum plane and must disap- 
pear in the full expression for the amplitude at each fixed order of approxi- 
mation. Using the result for the meson NLA impact factor such cancellation 
was demonstrated explicitly in Ref. [15] for the process in question. 



3. The impact factor for the virtual photon to light vector meson 

transition 

The definition of impact factor (IF) has been given in Ref. [6] ; in the case 
of scattering of the particle A off a Reggeized gluon with momentum q\ , for 
transverse momentum A and singlet color representation in the t-channel, 
the IF has the form [9] 



$A-A'(<fi, A, s ) 



£o 



so 



E / ^(*A " ^ d Pf T A{f} ( r l< 



{/}, 



d D - 2 k 



—$fZ n A'(k, A, s )ICr rn (k, qi, A) In 



Born i 



k 2 (k — A) 



(3.1) 



Here co(t) is the Reggeized gluon trajectory in the LLA. The integration in 
the first term of Eq. (|3.ip is done over the phase space dpf and over the 
squared invariant mass k of the system {/} produced in the fragmentation 
region of the particle A, r^rya are the related particle- Reggeon effective 



vertices. The second term in Eq. (J3JJ) is the counterterm for the LLA part 
of the first one, so that the logarithmic dependence of both terms on the 
intermediate parameter sa — ► oo disappears in their sum. The scale so 
is artificial and must disappear in the amplitude, to the given accuracy. 
The definition (|3.ip guarantees the infrared finiteness of the IFs of colorless 
particles 



G 




Here we study the NLA forward (A = 0) IF for the transition of a 
virtual photon to a light neutral meson & y *^y, V = p°,lo,4> (see Fig. 1). 
We use the auxiliary Sudakov vectors p\ and P2, such that p\ = p\ = and 
2(piP2) = s. The virtual photon momentum is p = pi — (Q 2 /s)p2, while 
Reggeon momenta are 

k + Q 2 + q 2 2 2-2 

Q = P2 + q±, q = q± = -q , 

s 

q' = i ^~V2 + q± ■ (3.2) 

s 

In the forward case under consideration the momentum transfer vector has 
only the longitudinal component, q — q' = CP2- Both the square of the 
Reggeon transverse momentum q 2 and the virtuality of the photon Q 2 are 
assumed to be much larger than any hadronic scale. Thus we neglect all 
power suppressed contributions and therefore the mass of the vector meson 
mass is put equal to zero and its momentum is identified with the Sudakov 
vector p\. 

It is possible to show (see Ref. |15| for details) that in this kinematics 
the IF can be calculated in the collinear factorization framework |25|, [26l [27] 
which was developed for the QCD description of hard exclusive processes. 
It turns out that the dominant helicity amplitude is a transition of the 
longitudinally polarized photon into the longitudinally polarized meson 
Vli and that both in LLA and in NLA the expression for the IF factorizes 
into the convolution^ 

4 5 CC ' 1 

^r L ^v L (a, s ) = — %& ^f V Q — / dzT H (z,a,s ,HF,^R)<P\\(z,HF) (3.3) 



2 Here and in the following we consider the color unprotected IF. 



7 



of a perturbatively calculable hard-scattering amplitude, Th, and a meson 
twist- 2 distribution amplitude, 4>\\(z) [28]. Here fjL 2 F ~ Q 2 ,q 2 is a factoriza- 
tion scale at which soft and hard physics factorizes according to Eq. (|3.3|) . 
The variable z corresponds to the longitudinal momentum fraction carried 
by the quark, for the antiquark the fraction is z = 1 — z. Finally, we intro- 
duced the ratio a = q 2 /Q 2 - For the cases of p°, u> and <j> meson production, 
e q in Eq. (|3,3p should be replaced by e/v2, e/(3\/2) and — e/3, respectively. 

We perform calculations with unrenormalized quantities, the bare strong 
coupling constant and the bare meson distribution amplitude 0y (z). 
Therefore the NLA expression for the hard-scattering amplitude Th ex- 
pressed in terms of these quantities contains both ultraviolet and infrared 
divergences, appearing as poles in the common dimensional regularization 
parameter e. The ultraviolet divergences disappear after the strong coupling 
constant renormalization. The surviving infrared divergences are only due 
to collinear singularities, the soft singularities cancel as usual after summing 
the "virtual" and the "real" parts of the radiative corrections. Since IFs 
should be infrared-finite objects for physical transitions, it must be possible 
to absorb the remaining infrared divergences into the definition of the non- 
perturbative distribution amplitude. This is achieved by the substitution of 
the bare distribution amplitude by the renormalized one (see Ref. [TS] for 
details) . 

The calculation of the lowest order contribution to the IF is straightfor- 
ward (see Ref. [15] ) . For leading power asymptotics, the dominant contri- 
bution is given by the production of a longitudinally polarized meson. The 
LLA result for the hard-scattering amplitude entering the IF is 

TP(z,a,s ,fi F ,fi R ) = a s — ^-^ . (3.4) 

a + zz 

Due to the collinear factorization, which effectively puts some fermion lines 
on the mass-shell, the complexity of the intermediate state contributing to 
the IF is reduced in comparison to the case of the virtual photon IF <I> 7 *^ 7 *. 
Actually things go as if we have one particle less in the intermediate state. 

In the NLA there are two contributions to the IF, from the two particle 
quark-antiquark (qq) and from three-particle quark-antiquark-gluon (qqg) 
intermediate states: 

T Q) = T (qq) + T (qqg) _ (3.5) 

To calculate the IF in the NLA one has to know the (qq) production vertices 
with NLA accuracy and the (qqg) ones at the Born level. To calculate 

one needs to convolute the NLA photon-Reggeon vertex ^^* q q [E] with the 
Born Reggeon-meson vertex Yyl - [T5] and the Born photon-Reggeon vertex 
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[291 03] with the NLA Reggeon-meson vertex T { yl q - [15]. To evaluate 
j^img) one needs to convolute Born photon-Reggeon vertex Ty?^ [13] with 

the Born Reggeon-meson vertex Fylqqg P^l - 

Taking also into account the so-called SA-counterterm (last term in 
Eq. (|3.ip ) and summing up all the contributions, one gets for the renor- 
malized hard scattering amplitude 

T H (z,a, s , Hf, m)\a-^o = us(hr) — (l + as ^ R ^ [t(z) + t(z)] \ , 

a + zz { 47T J 

(3.6) 

where the expression for r(z) is given in Ref. [15j . 

Using this result and the forward BFKL Green's function it is possible 
to build the forward amplitude of the process 7*(Qi)72(Q|) — ^ P1P2 i n the 
NLA. 



4. Building up the amplitude 

It is convenient to work in the transverse momentum representation, 
where "transverse" is related to the plane orthogonal to the vector mesons 
momenta. In this representation, defined by 

q Wi) = qM) , (4.1) 

(gi|<f 2 > = <5 (2) (<fi - q 2 ) , (A\B) = (A\k)(k\B) = J d 2 kA(k)B(k) , (4.2) 
the kernel of the operator K is 

K(q 2 ,q 1 ) = (q 2 \K\q 1 ) (4.3) 
and the equation for the Green's function reads 

l = (u- K)^ , (4.4) 

its solution being 

6 u = (u- K)- 1 . (4.5) 
The kernel is given as an expansion in the strong coupling, 

K = a s K° + alk 1 , (4.6) 

where 

a s = (4.7) 

7T 
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and N c is the number of colors. In Eq. (|4.6j) K° is the BFKL kernel in the 
LLA, K 1 represents the NLA correction. 

The impact factors are also presented as an expansion in a s 

where fy is the meson dimensional coupling constant (f p ~ 200 MeV) and 
e q should be replaced by e/v2, e/(3\/2) and — e/3 for the case of p°, lo and 
</> meson production, respectively. 

In the collinear factorization approach the meson transition impact fac- 
tor is given as a convolution of the hard scattering amplitude for the pro- 
duction of a collinear quark-antiquark pair with the meson distribution 
amplitude (DA). The integration variable in this convolution is the fraction 
z of the meson momentum carried by the quark (z = 1 — z is the momentum 
fraction carried by the antiquark): 







The NLA correction to the hard scattering amplitude, for a photon with 
virtuality equal to Q 2 , is defined as follows 

C«(<f 2 ) = ^ / dz f [t(z) + r(l - z)] h (z) , (4.10) 
c o 

with t{z) given in the Eq. (75) of Ref. [15]. C^\(q 2 ) are given by the 
previous expression with Q 2 replaced everywhere in the integrand by Q\ 
and Q?,, respectively. 

The distribution amplitude may be presented as an expansion in Gegen- 
bauer polynomials 

4>\\{z^ F ) = Qz{l-z) [l + a 2 {^ F )Cl l2 {2z-l) + a A {fi F )Cl /2 {2z-l) + ... . 

(4.11) 

The scale dependence of a n (/ip) is well known [251 l26j |2~7]: 

a n (p F ) = L^/P Q a n (fi) , (4.12) 
where L = a s (fi F )/a s (n) and 

lliV c 2n f 
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is the leading coefficient of the QCD /3-function, with nj the number of 
active quark flavors. The anomalous dimensions 7„ are positive and grow 
with n. Therefore any DA approaches the asymptotic form (p^(z) = 6z(l — 

z) at large 

Below we will use the DA in the asymptotic form. Besides the simplicity 
of the following presentation, the reason is twofold. Presumably, the form 
of DA chosen at low p p will affect mainly only the overall normalization of 
the amplitude but not the sum of BFKL energy logarithms and the result- 
ing dependence of the amplitude on the energy in which we are primarily 
interested in this study. Another point is that, according to QCD sum rules 
estimates [2H], a 2 (l GeV) is 0.18 ± 0.10 for p and ± 0.1 for 0. Therefore 
cj)^ may be indeed a good approximation for the DA of light vector mesons. 
Integrating over z in (j4.9j) with <j)\\(z, pyp) = (j)^(z), we obtain, for photon 
virtuality Q 2 , 



(a = C) = 6a 



Q' 2 



a . 2c + 1 
1 — — m 

c 2c- 1 



(4.14) 



where c = \fot + 1/4 . c[°l are given by the previous expression with Q 2 

replaced by Q 2 and Q?,, respectively. For the NLA term C^liq 2 ) the inte- 
gration over z can be performed by a numerical calculation. 

To determine the amplitude with NLA accuracy we need an approximate 
solution of Eq. (|4.5p . With the required accuracy this solution is 



G u = (u- ask )- 1 + {uj- ask )' 1 (alk 1 ^ {uj - a^ )" 1 + O (o^K 1 

(4.15) 

The basis of eigenfunctions of the LLA kernel, 
K°\u) = x (v)\u) , X iy) = 2^(1) - i> Q + ii^j - ill Q - if) , (4.16) 

is given by the following set of functions: 

^ = jti (^r 4 • <4 - i7) 

for which the orthonormality condition takes the form 

= / H ( g ) =5(u-S). (4.18) 



3 The dependence of the resulting amplitude on is subleading. Due to the collinear 
counterterm, see Eq. (72) of [15] . the NLA correction to the meson impact factor 
contains a term proportional to \n(fiF), see Eq. (75) of [15], which compensates in 
the amplitude with NLA accuracy the effect of the meson DA variation with /if. 
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The action of the full NLA BFKL kernel on these functions may be expressed 
as follows: 



A'!") - <Ui>h)\<^\") + «:'(///,>) ( x (1 V) + ^XMM^) ) I'') 



00 



(4.19) 



where the first term represents the action of LLA kernel, while the second 
and the third ones stand for the diagonal and the non-diagonal parts of the 
NLA kernel. The function x { v \ calculated in Ref. [2], is conveniently 
represented in the form 



x (1) H 



A) 

8N, 



(4.20) 



where 



XW 



7T 



2 -4 



- 6C(3) - X "{v) 



7T 



COsh(7TZ/) 



7T 2 sinh(^) / / n£\ \\ + Ylv l \ 
+ 2z,cosh 2 (™) I + V m) 16(1 + ^ 2 ) J + 



(4.21) 



4>(y) = 2 dx 



cos(^ ln(x)) 
[l+x)^x~ 



7T 



Li9(xj 

6 v ; 



Li 2 (x) 



dt 



ln(l - i) 



(4.22) 



Here and below x'( u ) = d(x(v))/dv and x"( v ) = ^ 2 (x( u )) I d 2 v . 

We will need also the \v) representation for the impact factors, which is 
defined by the following expressions 



c W(a 2 ) + r° C {0) (a 2 ) + r° 

1 W j " ' dv' Cl {v>)(v'\<t) , ^-Jfi= / di/c 2 (^( g », 



(4.23) 



(0)^-2, (i!) 



•2\f-3 



J IT 



7T 



V2 ' 
(4.24) 

and by similar equations for c^(y) and (^'(y) from the NLA corrections 
to the impact factors, (q 2 ) and C^\q 2 ). 
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Using (I4.15P and (I4.19P one can derive, after some algebra, the following 
representation for the amplitude 



Xm s (A) s 
B X D 2 = (2vr) 2 



+00 



dv 



•so 



a«(/iji)xM 



o? s (hr)c 1 (v)c 2 {v) 



1 + a s (ji R ) 



+ 8]v: x( " ) 



r, ( /-') Co ( /' 



Cl(l/ 



c 2 (z/) 
, , 10 . 



+ Q^(/x R )ln 



x(f) (4-25) 



10 dH'^: , ol 2 



21n(/i 



We find that 



c l,2^) 



Ql,2 



zr.iv- 



r 2 [f±ii/] 6vr 



c 1 (i/)c 2 (i/) 



r[3 ± 2iv] cosh(vr^) ' 
1 (QlY 9vr 3 (l + 4i/ 2 )sinh(7ri/) 



Q1Q2 \QiJ 32 v (1 + z/ 2 ) cosh^TTz/) 



(4.26) 
(4.27) 



of// 



V>(3 + 2zV) +^(3 - liv) 



(4.28) 



^ I — h iv — ip 



if ) - In (Q1Q2 



It can be useful to separate from the NLA correction to the impact factor 
the terms containing the dependence on sq and on /3q, 



C^{q 2 ) = / dz 



q 2 + zzQ 2 



0||(*) 



(4.29) 



hi 



so 

Q 2 



In 



(a + zz) 4 
a 2 z 2 z 2 



+ 



In 



+ 3 - ln(a) j + . . . 



Accordingly, one can write 

4>) = ci>)+4>) 



(4.30) 
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where c{\(v) are the contributions from the terms isolated in the previous 

equation and c^\{v) represent the rest. After straightforward calculations 
we found that 



d{v) c 2 (v) 



Q1Q2 



2N r 



In 



/4 



5 

Q1Q2J + 3 



(4.31) 



+ V(3 + 2iu) + ij>(3 - 2iu) -il>\^ + iv)-if>0; 
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Using Eq. (|4.25p we construct the following representation for the am- 
plitude 



QiQ 2 Tm s A 1 
D X D 2 s ~ (2^) 



2 a s (ji R ) 2 



b + V a s (n R ) n b n ( In ( — ) + d n (s , fj, R ) In ( — 
^ V \s J \s 



(4.32) 



where the coefficients 



+00 



Q1Q2 



dv c\{y)c%{y) 



nl 



(4.33) 



are determined by the kernel and the impact factors in LLA. Note that 



bo = ^ (7C(3) - 6) , 



(4.34) 



therefore in the Born (the 2-gluon exchange) limit our result coincides with 
that of Ref. [18]. 
The coefficients 



d n = n In 



-so 



n(n — 1) 



4A, 



Q1Q2 



+ 



Q1Q2 



dv (n + l)f(v)c 1 (v)c 2 (v) 



X n ">) 
(n- 1)! 



(4.35) 



+ 



Q1Q2 



c\{v)c 2 {v) 



X n ">) 
(n- 1)! 



cfV) , 4 1 V) | (n 



ciH c 2 {v) 
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are determined by the NLA corrections to the kernel and to the impact 
factors. Here we use the notation 



f{y) = - + ^(3 + 2w) + if>(3 - 2iv) - V (j + iA - il> ( 



IV 



(4.36) 



One should stress that both representations of the amplitude (I4.32P 
and (|4.25p are equivalent with NLA accuracy, since they differ only by 
next-to-NLA (NNLA) terms. Actually there exist infinitely many possi- 
bilities to write a NLA amplitude. For instance, another possibility could 
be to exponentiate the bulk of the kernel NLA corrections 



Im s (A) 



(2,? J dV U 

— oo 



x Oi 2 s (n R )ci{v)c 2 {v) 



1 + a s (n R 



a 2 s ((j, R ) In — 



so J 8N C 




(4.37) 



This form of the NLA amplitude was used in Ref. [30] (see also [31J), without 
account of the last two terms in the second line of (|4.37p . for the analysis 
of the total 7*7* cross section. 

Since as we will shortly see the NLA corrections are very large, the choice 
of the representation for the NLA amplitude becomes practically important. 
In the present situation, when an approach to the calculation of the NNLA 
corrections is not developed yet, the series representation (|4.32|) is, in our 
opinion, a natural choice. It includes in some sense the minimal amount of 
NNLA contributions; moreover, its form is the closest one to the initial goal 
of the BFKL approach, i.e. to sum selected contributions in the perturbative 
series. 

It is easily seen from Eqs. (|4.32p - (|4.36p that the amplitude is independent 
in the NLA from the choice of energy and strong coupling scales. Indeed, 
with the required accuracy, 




(4.38) 



and therefore terms a"ln n_1 slnso and a^ln" --1 sln///j cancel in (I4.32p . 

One can trace the contributions to each d n coefficient coming from the 
NLA corrections to the BFKL kernel and from the NLA impact factors 



(4.39) 
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jIF 



nln 



so 




(4.40) 



The first coefficient, d\, is entirely due to the NLA corrections to the impact 
factors, 



d x = < , d\ er = . (4.41) 

Let us note that in the BFKL formalism the NLA contribution to the impact 
factors guarantees not only independence of the amplitude from the energy 
scale, so, but it contains also a term proportional to In hr which is important 
for the renorm-invariance of the predicted results, i.e. the dependence of 
the amplitude on fiR and sq is subleading to the NLA accuracy. 



5. Numerical results 

In this Section we present some numerical results for the amplitude 
given in Eq. (|4.32p for the Qi = Q2 = Q kinematics, i.e. in the "pure" 
BFKL regime. The other interesting regime, Qi 3> Q2 or vice-versa, where 
collinear effects could come heavily into the game, will not be considered 
here. We will emphasize in particular the dependence on the renormalization 
scale fj,R and so in the NLA result. 

In all the forthcoming figures the quantity on the vertical axis is the 
L.H.S. of Eq. (|4.32p . Zm s (A)Q 2 / (s D1D2). In the numerical analysis pre- 
sented below we truncate the series in the R.H.S. of Eq. (|4.32|) to n = 20, 
after having verified that this procedure gives a very good approximation 
of the infinite sum for the Y values Y < 10. We use the two-loop running 
coupling corresponding to the value a s {Mz) = 0.12. 

We have calculated numerically the b n and d n coefficients for nj = 5 
and s = Q 2 = n 2 R , getting 

60 = 17.0664 

61 = 34.5920 6 2 = 40.7609 63 = 33.0618 64 = 20.7467 
65 = 10.5698 b 6 = 4.54792 b 7 = 1.69128 6 8 = 0.554475 



dt = -3.71087 d 2 = -11.3057 d 3 = -23.3879 d A = -39.1123 
d 5 = -59.207 d 6 = -83.0365 d 7 = -111.151 d 8 = -143.06 . 

(5.1) 
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In this case contributions to the d n coefficients originating from the NLA 
corrections to the impact factors are 

df = -3.71087 df = -8.4361 4 F = -13.1984 df = -18.0971 
df = -23.0235 df = -27.9877 df = -32.9676 df = -37.9618 . 

(5.2) 

Thus, comparing (15. ip and (|5.2p . we see that the contribution from the 
kernel starts to be larger than the impact factor one only for n > 4. 

These numbers make visible the effect of the NLA corrections: the d n 
coefficients are negative and increasingly large in absolute values as the 
perturbative order increases. The NLA corrections turn to be very large. 
In this situation the optimization of perturbative expansion, in our case the 
choice of the renormalization scale p,R and of the energy scale sq, becomes 
an important issue. Below we will adopt the principle of minimal sensitivity 
(PMS) |32j . Usually PMS is used to fix the value of the renormalization 
scale for the strong coupling. We suggest to use this principle in a broader 
sense, requiring in our case the minimal sensitivity of the predictions to the 
change of both the renormalization and the energy scales, p,R and sq. 

Since the dependence of results on sq is a feature typical of the BFKL 
approach and is somewhat new for the application of PMS, we will first 
illustrate the success of PMS in this respect on the following QED result 
known since a long time. In 1937 Racah calculated the total cross section 
for the production of e + e~ pairs in the collisions of two heavy ions at high 
energies [33j, 

-%f% + ^ + c)+o(^); (5.3, 

here Zip are the ions charges, m e is electron mass, the ions' four-momenta 
are p l>2 , 

Z = ln^iM, (5.4) 
mim% 

is the energy logarithm and mi 2 are the masses of the ions. The contribu- 
tions suppressed by the power of energy are denoted as 0(l/(pip2))- 

The coefficients in front of the subleading logarithms are large and have 
alternating signs 

A = -178/28 = -6.35714 

B = -"-(77T 2 + 370) = 15.6817 (5.5) 
28 

1/13 \ 
C = -— (348 + yvr 2 - 21 C (3) = -13.8182 . 
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O/Oq 




Fig. 3. cr/cro as a function of the energy logarithm I for the cases of exact result 
of Racah, approximated result with l = —A/3 (PMS optimal choice) and Iq = 
(kinematical scale for energy logarithms). 



To illustrate the application of PMS, imagine that we know only the coeffi- 
cient A in front of the first subleading logarithm. Then using this knowledge 
we can construct the following approximation 

a^ = a ((l-l f + (A + 3l )(l-lo) 2 ) , *o = ^^J^ ■ ( 5 - 6 ) 

(an analog of NLA in the BFKL approach) where we shift the energy scale 
introducing the parameter Iq. Note that the dependence of the cross section 
on Iq, is subleading in the approximation used in Eq. (|5.6p . We fix Iq by 
requiring the minimal sensitivity of (|5,6p to the change of this parameter. 
It is not difficult to find that this procedure gives l = -A/3 = 2.119050 In 
Fig. Owe present three curves for a /gq as a function of the energy logarithm 
I; the first one was calculated using the exact result of Racah (with all 
subleading logarithms), the other two curves were calculated using (|5.6p 
with l = and with the PMS value Iq = 2.11905. 

We see that the PMS approach gives a very good approximation to the 
Racah resul10. On the other hand the procedure with Iq = 0, which means 
that a kinematical scale for energy logarithms is used in the approximate 
formula, makes an awfully bad job for the whole I range presented in the 
figure. 

4 Note that in this example PMS gives the value of the parameter lo for which the 
correction to the lowest approximation, (J— lo) s , vanishes. Therefore in this case PMS 
gives a result which coincides with the one given by another alternative approach to 
optimize the approximation, the fast apparent convergence prescription |34] , 
The negative cross section at I < 2 is due to the fact that terms subleading in energy, 
0(l/(p\p2)) in (|5.3[) . are not taken into account. 
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Fig. 4. (Left) lm s (A)Q 2 /(sD 1 D 2 ) as a function of Y a at //r = 10Q. The different 
curves are for Y values of 10, 8, 6, 4 and 3. The photon virtuality Q 2 has been 
fixed to 24 GeV 2 (n/ = 5). (Right) lm s (A)Q 2 / (s DiD 2 ) as a function of at 
Y=6. The different curves are, from above to below, for Yq values of 3, 2, 1 and 0. 
The photon virtuality Q 2 has been fixed to 24 GeV 2 (nf = 5). 

Returning to our problem, we apply PMS to our case requiring the min- 
imal sensitivity of the amplitude (|4.32p to the variation of fiR and Sq. More 
precisely, we replace in (|4,32p ln(s/so) with Y — Yq, where Y = ln(s/Q 2 ) 
and Yq = \~h(sq/Q 2 ), and study the dependence of the amplitude on Yq. 

The next two figures illustrate the dependence on these parameters for 
Q 2 =24 GeV 2 and nj = 5. In Fig. [Deleft) we show the dependence of ampli- 
tude on Yq for /ir = WQ, when Y takes the values 10, 8, 6, 4, 3. 

We see that for each Y the amplitude has an extremum in Yq near which 
it is not sensitive to the variation of Yq, or sq. Our choice of jir for this 
figure is motivated by the study of {j,r dependence. In Fig. fright) we 
present the hr dependence for Y = 6; the curves from above to below are 
for y =3, 2, 1, 0. 

Varying fiR and Yq we found for each Y quite large regions in jiR and Yq 
where the amplitude is practically independent on fiR and Yq. We use this 
value as the NLA result for the amplitude at given Y. In Fig. [5] we present 
the amplitude found in this way as a function of Y. The resulting curve is 
compared with the curve obtained from the LLA prediction when the scales 
are chosen as /j,r = 10Q and Yq = 2.2, in order to make the LLA curve 
the closest possible (of course it is not an exact statement) to the NLA one 
in the given interval of Y. The two horizontal lines in Fig. [5] are the Born 
(2-gluon exchange) predictions calculated for [Ir = Q and \xr = 10Q. 

Similar procedure was applied to a lower value of the photon virtuality, 
Q 2 =5 GeV 2 and n f = 4 (see Ref. [E] for details). 

We stress that one should take with care BFKL predictions for small 
values of Y, since in this region the contributions suppressed by powers of 
the energy should be taken into account. At the lowest order in a s such con- 
tributions are given by diagrams with quark exchange in the t-channel and 
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Fig. 5. lm s (A)Q 2 / (s DiD 2 ) as a function of Y for optimal choice of the energy 
parameters Y and fiR (curve labeled by "NLA"). The other curves represent the 
LLA result for Y — 2.2 and //r = 10Q and the Born (2-gluon exchange) limit 
for \xr = Q and ixr — 10Q. The photon virtuality Q 2 has been fixed to 24 GeV 2 
(n f = 5). 



are proportional in our case to aEM&sfv/Q 2 - At higher orders power sup- 
pressed contributions contain double logarithms, terms ~ a™ln 2n s, which 
can lead to a significant enhancement. Such contributions were recently 
studied for the total cross section of 7*7* interactions [35]. 

If the NLA (and LLA) curves in Fig. [5] are compared with the Born 
(2-gluon exchange) results, one can conclude that the summation of BFKL 
series gives negative contribution to the Born result for Y < 6 if one chooses 
for the scale of the strong coupling in the Born amplitude the value given 
by the kinematics, = Q. We believe that our calculations show that 
one should at least accept with some caution the results obtained in the 
Born approximation, since they do not give necessarily an estimate of the 
observable from below. 

Another important lesson from our calculation is the very large scale for 
a s (and therefore the small a s itself) we obtain using PMS. It appears to 
be much bigger than the kinematical scale and looks unnatural since there 
is no other scale for transverse momenta in the problem at question except 
Q. Moreover one can guess that at higher orders the typical transverse 
momenta are even smaller than Q since they "are shared" in the many-loop 
integrals and the strong coupling grows in the infrared. In our opinion the 
large values of fiR we found is not an indication of the appearance of a new 
scale, but is rather a manifestation of the nature of the BFKL series. The 
fact is that NLA corrections are large and then, necessarily, since the exact 
amplitude should be renorm- and energy scale invariant, the NNLA terms 
should be large and of the opposite sign with respect to the NLA. We guess 
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that if the NNLA corrections were known and we would apply PMS to the 
amplitude constructed as LLA + NLA-corrections + NNLA-corrections, we 
would obtain in such calculation more natural values of fiR. 

In the last years strong efforts have been devoted to the improvement of 
the NLA BFKL kernel as a consequence of the analysis of collinear singu- 
larities of the NLA corrections and by the account of further collinear terms 
beyond NLA [21, 36J. This strategy has something in common with ours, 
in the sense that it is also inspired by renormalization-group invariance and 
it also leads to the addition of terms beyond the NLA. These extra-terms 
are large and of opposite sign with respect to the NLA contribution, so that 
they partially compensate the NLA corrections. The findings of the present 
work suggest, however, that the corrections to the impact factors heavily 
contribute to the NLA amplitude, being even dominating in some interval of 
non-asymptotically high energies. Moreover, by inspection of the structure 
of the amplitude in the regime of strongly asymmetric photon virtualities, 
one can deduce that also the impact factors generate collinear terms which 
add up to those arising from the kernel, see e.g. Eqs. (84) and (85) of 
Ref. [15] . This leads us to the conclusion that in the approaches based on 
kernel improvement the additional information coming from impact factors 
should somehow be taken into account when available. These issues cer- 
tainly deserve further investigation and we believe that useful hints in this 
direction can be gained from the study of the 7*7* — ► VV amplitude in the 
regime of strongly ordered photon virtualities )22[ [23] . 

We conclude this Section with a comment on the possible implications 
of our results for mesons electroproduction to the phenomenologically more 
important case of the 7*7* total cross section. By numerical inspection 
we have found that the ratios b n /bo we got for the meson case agree for 
n = 1 -J- 10 at 1 -T- 2% accuracy level with the analogous ratios for the 
longitudinal photon case and at 3.5 4- 30% accuracy level with those for 
the transverse photon case. Should this similar behavior persist also in the 
NLA, our predictions could be easily translated to estimates of the 7*7* 
total cross section. 

6. Study of systematic effects 

It is important to have an estimate of the systematic uncertainty which 
plagues our determination of the energy behavior of the amplitude. The 
main sources of systematic effects are given by the choice of the represen- 
tation of the amplitude and by the optimization method adopted. In the 
following, we compare the determination of the amplitude at Q 2 = 24 GeV 2 
(rif = 5) through the PMS method, given in Fig. [5l with other determina- 
tions obtained changing either the representation of the amplitude or the 
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Fig. 6. (Left) Xm s {A)Q 2 /{sD 1 D 2 ) as a function of Y at Q 2 =24 GeV 2 (n f = 5) 
from series and "exponentiated" representations, in both cases with the PMS opti- 
mization method. (Right) Xm s (A S eries)Q 2 / (s D1D2) as a function of Y at Q 2 =24 
GeV 2 (nj = 5) from the series representation with PMS and FAC optimization 
methods. 



optimization method. 

At first, we compare the series and the "exponentiated" determinations 
using in both case the PMS method. The procedure we followed to deter- 
mine the energy behavior of the "exponentiated" amplitude is straightfor- 
ward: for each fixed value of Y we determined the optimal choice of the 
parameters hr and Yq for which the amplitude given in Eq. (|4.37p is the 
least sensitive to their variation. Also in this case we could see wide regions 
of stability of the amplitude in the (fJ,R, Yq) plane. The optimal values of fiR 
and Yq are quite similar to those obtained in the case of the series represen- 
tation, with only a slight decrease of the optimal hr. In Fig.^left) we show 
the result and compare it to the PMS determination from the series repre- 
sentation. The two curves are in good agreement at the lower energies, the 
deviation increasing for large values of Y. It should be stressed, however, 
that the applicability domain of the BFKL approach is determined by the 
condition a s (fj,R)Y ~ 1 and, for Q 2 =24 GeV 2 and for the typical optimal 
values of hr, one gets from this condition Y ~ 5. Around this value the 
discrepancy between the two determinations is within a few percent. 

As a second check, we changed the optimization method and applied it 
both to the series and to the "exponentiated" representation. The method 
considered is the fast apparent convergence (FAC) method [33], whose strat- 
egy, when applied to a usual perturbative expansion, is to fix the renormal- 
ization scale to the value for which the highest order correction term is 
exactly zero. In our case, the application of the FAC method requires an 
adaptation, for two reasons: the first is that we have two energy parameters 
in the game, fiR and Yq, the second is that, if only strict NLA corrections 
are taken, the amplitude does not depend at all on these parameters. 

Therefore, in the case of the series representation, Eq. (|4.32p . we choose 
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Fig. 7. (Left) 2m s (A)Q 2 /(sD 1 D 2 ) as a function of Y at Q 2 =24 GcV 2 (n f = 
5) from the series representation with the PMS optimization method and from 
the "exponentiated" representation with the FAC optimization method. (Right) 
lm s (A)Q 2 / (s DiD 2 ) as a function of Y at Q 2 —24 GeV 2 (rif = 5) from the series 
representation with PMS and BLM optimization methods. 



to put to zero the sum 
1 



( 27r ) £i 



n-l 



and found for each fixed Y the values of hr and lo f° r which the vanishing 
occurs. This gives a line of values in the (fJ,R, Yq) plane, among which the 
optimal choice is done applying a minimum sensitivity criterion. The result 
is shown in Fig. fright). The agreement with the series representation with 
the PMS method is rather good over a wide energy range. 

In the case of the "exponentiated" amplitude" , Eq. (|4.37[) , we proceeded 
in the same way, but requiring the vanishing of the expression given by the 
R.H.S. of Eq. (|Q7|) minus the LLA amplitude, i.e. 



Im s (A, 



exp J 



D 1 D 2 



(2nf 



dv 



c? s {iiR)cx[y)ci[y) . 



In Fig. [?K left) the result is compared with series representation in the PMS 
method: there is nice agreement over the whole energy range considered. 

Another popular optimization procedure is the Brodsky-Lepage-Mac- 
kenzie (BLM) method [37], which amounts to perform a finite renormaliza- 
tion to a physical scheme and then choose the renormalization scale in order 
to remove the /^-dependent part. We applied this method only to the series 
representation, Eq. (|4.32|) . and proceeded as follows: we first performed a 
finite renormalization to the momentum (MOM) scheme with £ = (see 
Ref. 
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Fig. 8. (Left) Linear-log plot of da/dt\ t=t [pb/GeV 2 ] as a function of y/s at Q 2 = 16 
GeV 2 (n / = 4) from the series representation with the PMS optimization method 
(solid line) compared with the determination from the approach in Ref. [19] (dashed 
line). (Right) The same as (left) at Q 2 =4 GeV 2 (n f = 3) 
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then, we chose Yq and pr in order to make the term proportional to (3q 
in the resulting amplitude vanish. We observe that the /^-dependence in 
the series representation of the amplitude is hidden into the d n coefficients, 
Eq. (|4.35p . Among the resulting pairs of values for Yq and pr, we deter- 
mined the optimal one according to minimum sensitivity. This method has 
a drawback in our case, since for each fixed Y, the optimal choice for Yq 
turned to be always Yq ~ Y. However, if one blindly applies the proce- 
dure above, one gets a curve which slightly overshoots the one for the series 
representation with the PMS method, see Fig. Upright). 



7. The differential cross section at the minimum \t\: comparison 
with an approach based on collinear improvement 

The 7*7* — > pp process at the lowest order (two-gluon exchange in the 
i-channel) was studied in Ref. [18]. At that level our results coincide, see 
also |16j . The same process with the inclusion of NLA BFKL effects has 
been considered in Ref. [TU]. In that paper, the amplitude has been built 
with the following ingredients: leading-order impact factors for the 7* — > p 
transition, BLM scale fixing for the running of the coupling in the prefactor 
of the amplitude (the BLM scale is found using the NLA 7* — > p impact 
factor calculated in Ref. [15J) and renormalization-group-resummed BFKL 
kernel, with resummation performed on the LLA BFKL kernel at fixed 
coupling [38] . In Ref. [19] the behavior of da/dt at t = to was determined 
as a function of y/s for three values of the common photon virtuality, Q=2, 
3 and 4 GeV. 
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Fig. 9. (Left) Linear plot of da / dt\ t=to [pb/GeV 2 ] as a function of y/s at Q 2 = 16 
GeV 2 (rif =4) from the series representation with the PMS optimization method 
using NLO impact factors (solid line) and LO impact factors (dashed line). (Right) 
The same as (left) at Q 2 =4 GeV 2 (n f = 3). 



In order to make a comparison with the findings of Ref. [19], we com- 
puted da/dt at t = to for Q=2 and Q=A GeV as functions of i/i. We 
used / p =216 MeV, oem = 1/137 and the two-loop running strong coupling 
corresponding to the value a s (Mz) = 0.12. The results are shown in the 
linear-log plots of Figs. [51 which show disagreement. This is not surprising 
in consideration of the approximations adopted in Ref. [19j , 

It would be interesting to understand to what extent this disagreement 
is due to the use in Ref. [19] of LLA impact factors instead of the NLA ones 
or to the way the collinear improvement of the kernel is performed. 

In order to understand to what extent the discrepancy is due to the use 
of leading order (LO) impact factors instead of next-to-leading order (NLO) 
ones, we repeated our determination of da/dt at t = to for Q=2 and Q=4 
GeV, using LO impact factors and keeping from the their NLO contribution 
only the terms proportional to lnfso/CQiQa)] and to ln[//|j/ (Q1Q2)] which 
are universal and needed to guarantee the so- and //^-independence of the 
amplitude with NLA accuracy. The result is that da/dt at t = to increases 
roughly by an order of magnitude with respect to our previous determination 
(see Figs. [9|) and therefore the disagreement with [19] becomes even worse. 
This is not surprising: impact factors give a sizable contribution to the NLA 
part of the amplitude which is negative with respect to the LLA part; if they 
are kept at LO, the NLA part of the amplitude is less negative and the total 
amplitude is therefore increased. 



8. Conclusions 

We have determined the amplitude for the forward transition from two 
virtual photons to two light vector mesons in the Regge limit of QCD with 
next-to-leading order accuracy. This amplitude is the first one ever written 
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in the next-to-leading approximation for a collision process between strongly 
interacting colorless particles. It is given as an integral over the v parame- 
ter, which labels the eigenvalues of the leading order forward BFKL kernel 
in the singlet color representation. This form is suitable for numerical eval- 
uations. The result obtained is independent on the energy scale sq, and on 
the renormalization scale /j,r within the next-to-leading approximation. 

Using a series representation of the amplitude which includes the depen- 
dence on the energy scale and on the renormalization scale at subleading 
level, we performed a numerical analysis in the kinematics when the two 
colliding photons have the same virtuality, i.e. in the "pure" BFKL regime. 
We have found that the next-to-leading order corrections coming from the 
kernel and from the virtual photon to light vector meson impact factors are 
both large and of opposite sign with respect to the leading order contribu- 
tion. 

An optimization procedure, based on the principle of minimal sensitivity 
method, has proved to work nicely and has lead to stable results in the 
considered energy interval, which allows us to predict the energy behavior of 
the forward amplitude. The procedure consists in evaluating the amplitude 
at values of the energy parameters for which it is the least sensitive to 
variations of them. We have found that there are wide regions of values of 
so and [Ir where the amplitude remains almost flat. 

The optimal choices of sq and [ir are much larger that the kinemati- 
cal scales of the problem. More than being the indication of appearance 
of another scale in the problem, this could be related to the nature of the 
BFKL series. The renorm- and energy scale invariance, together with the 
large next-to-leading approximation corrections, call for large next-to-next- 
to-leading order corrections, which are most probably mimicked by unnat- 
ural optimal values for sq and ixr. 

The use of other optimization methods and/or different, equivalent in 
the NLA, representations of the amplitude gives results does not change the 
behaviour of the amplitude with energy and allows for an estimate of the 
systematic uncertainty of our determinations. 
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